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We show that Green function methods can be straightforwardly applied to nonlinear equations 
appearing as the leading order of a short time expansion. Higher order corrections can be then 
computed giving a satisfactory agreement with numerical results. The relevance of these results 
relies on the possibility of fully exploiting a gradient expansion in both classical and quantum field 
theory granting the existence of a strong coupling expansion. Having a Green function in this regime 
in quantum field theory amounts to obtain the corresponding spectrum of the theory. 
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I. INTRODUCTION 



Nonlinear equations represent a class of very difficult mathematical problems to manage by analytical methods. A 
lot of fundamental aspects of physics are described by these equations making not easy their understanding due to 
the lack of useful techniques. 

In this paper we present an approach that is based on an unexpected result for Green functions. First hints 

, , in this direction were obtained in [l], [2j where we showed that for a part of the integration interval, a nonlinear 

4h : differential equation like cj> + 4> z = j, being j a source term, can be solved by its Green function G + G 3 — S(t) as 

Although this solution was put forward, the knowledge of this result is hardly useful unless we are not able to 
understand how to get higher order corrections. The aim of this paper is to give a proper understanding of this 
i—i, solution and to give a technique to get higher order corrections in order to improve it. We will show that it represents 
a short time solution. Then, a form factor described by polynomial terms in time can correct properly the propagator 
\ to improve in some cases this approximation. 

The reason to give such a solution relies on the possibility to treat strong coupled quantum field theories that at the 
leading order produce nonlinear equations driven by a source. Strongly coupled theory can be managed by a gradient 
. expansion that is the dual perturbation series to a weak coupling expansion as we proved in [J H, [j] . By "dual" we 
t-H ' mean that two series can be obtained by simply interchanging the terms of the expansion producing in a case a series 
, with an expansion parameter being the inverse of the expansion parameter of the other series, the asymptotic series 
so obtained holding in the proper limit where this parameter gives a converging expansion (coupling going to infinity 
in a case while going to zero in the other). This approach is true for any differential equation set and we applied it 
also in general relativity [|| obtaining a sound proof of the Belinski-Khalatnikov-Lifshitz conjecture [1, @, 0] as this is 
| a result of a gradient expansion. 
• »"H , A strongly coupled system in quantum mechanics is known to be a classical system as was firstly shown by Simon 
Q. We revised this approach in Q where we have seen that the gradient expansion for the Schrodinger equation, 
also known as Wigner-Kirkwood expansion, gives rise to a Thomas-Fermi approximation to the leading order for 
a many-body system [101 ] and has the same eigenvalue expansion as for a WKB approximation. Wigner-Kirkwood 
expansion is indeed the gradient expansion of the Schrodinger equation. 

Gradient expansions in quantum field theories were not widely used before while their proper understanding is not 
that easy. Our aim in this paper is to fully exploit this perturbation approach and its application in quantum field 
theory wherever possible. This method may pave the way to manage analytically some problems that now appear 
difficult to manage also in a wide variety of fields where nonlinear equations are at the foundations. 

The paper is structured in the following way. In section |TT] we show how to derive a gradient expansion out of a 
duality principle in perturbation theory with the proper understanding of the expansion parameter. In section IIIII 
we present the main motivation for this paper, that is the continuum limit of a scalar quantum field theory giving 
rise to a model nonlinear equation we will use throughout the paper. In section IIVI we present the method firstly 
applied to a simple Riccati equation having a known analytical solution and then we generalize our method to the 
case of the leading order equation of a scalar field theory. In section [V] we give the numerical results showing how the 
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approximation improves with higher order corrections varying also the forcing term into the equation. In section lVTl we 
compare our approach with functional iteration method, a well known method used to analyze non-linear differential 
equations. This will give a proper understanding of the speed of convergence of our method. Finally, in section IVlTl 
conclusions are given. 



II. DUAL EXPANSION FOR NONLINEAR PDES 



In order to make the paper self-contained we present here some material already given in Ref. [i[ . We specialize the 
presentation to a A0 4 model that is our reference model. 
The Hamiltonian of the model is given by 



H= I d D ~ x x 



i7r 2 + i(V0) 2 + V(0) 



(1) 



being D the spacetime dimensionality and V((f>) = \<j) + j(f> and we take the case of a single component for the sake 
of simplicity. Hamilton equations are 

dt<t> = 7T (2) 
d t iT = V 2 0- 0- A0 3 . 

We can see at glance that we can chose to do perturbation theory by two different choices. One can take either A0 3 
or V 2 — as a small term. What we want to understand is the link between the two series with respect to the 
parameter A. 

By choosing A0 3 as a small term one gets the small perturbation series 

(3) 
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where it easily seen that the free theory, D0o + 0o = 0, is the leading order solution. Our aim is to derive a dual 
perturbation series to this one meaning by this that we want a series with a development parameter going as j . 
In order to reach our aim, following the principle of duality in perturbation theory we put 

T = V\t (4) 



7T = V\ 



1 1 

n ° + A 71 " 1 + A 2 n2 

^01 + ^2 + ... 



The following non trivial set of equations is obtained 

3 T (j) = 7T (5) 

<9 r 01 = TTl 
<9r0>2 = 7T2 



d T TTQ = -0o 

8 T TTl = V 2 00 - 00 - 30Q01 

d T ir 2 = V 2 0i - 0i - 30001 - 30Q02 
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whose solution proves the existence of a dual perturbation series for the classical Xcj) 4 theory. We easily realize that 
this set of equations would have been obtained if one takes as a small term V 2 </> — <f> giving rise in this case to a 
gradient expansion, that is a series having derivatives in space as small terms. So, strong coupling expansion and 
weak coupling expansion are related by the duality principle in perturbation theory llj producing in the former case 
a gradient expansion. This result can be easily generalized to any kind of PDE [1, H|. 

The point to be noted is that to have an analytical result for a strong coupling expansion we have to solve a 
nonlinear differential equation that in this case is given by 



d^o + ^ = 0. (6) 

Things can be more involved when a source term is present as is generally the case in quantum field theory and a 
meaning should be attached to the leading order equation 

d 2 M + $ = j. (7) 

being j a source term. The aim of this paper is to show that indeed an approach through Green functions is applicable 
in these cases, that is, as already shown by numerical methods in [J a first approximated solution is given by 

cj>= [ T dr'G(r-T')j(r') (8) 



JO 

being G(r) the Green function solving the equation 

d 2 G(r)+G(r)^S(r). (9) 

We will give in this paper a general approach to compute higher order corrections to this result. 
We will note that the method can be applicable when a solution is known to an equation like 

d 2 T G(r) + F(G(t)) = aS(r) (10) 

being a a constant and F(G(t)) a generic term. Otherwise we are not able to get analytical results and we have to 
resort to numerical methods. Anyhow, the situation is favorable for the most common models. 



III. GRADIENT EXPANSION AND QUANTUM FIELD THEORY 

Quantum field theory of the model we are considering is given by the generating functional 

Z\j] = J [^] e { l / dD -[5(a t 0) 2 -^0 4 +J0]} e {^/<i D ^^V0) 2 +i0 2 ]} (n) 

that we have written separating the spatial part from the rest. We did this in order to derive the strong coupling 
expansion to this case as already done in sec|TT] for the classical model. By doing the expansion, considering the 
gradient term i(V</>) 2 + as small, the leading order term to be computed is 

Z \j] = J [d0] e {*/ dI M^) 2 -^ 4 +^]} (12) 

and in the end we are left with the equation to solve 

dU + \<i?=3 (13) 

that is the leading order of our gradient expansion as already seen in secim 

The applicability of the Green function method implies that also in a strong coupling regime one can obtain 
information on the spectrum of the theory in this limit. We can exploit this point easily for a our case. Firstly, we 
use the mass /iq of the theory to make all adimensional putting x — > fiox, <fi 2 — > (J,q~ D ' <fi 2 and introducing the coupling 
constant g — . Then, let us consider the equation 

Mo 



d 2 G + gG 3 = 6(t) 



(14) 
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that has solution [l[ 



G(t) = 9(t) ( - 



sn 



(!)*<■' 



(15) 



being the Hcavisidc function and sn a Jacobi elliptical function. Being the equation second order we have that 
also the time reversed solution holds. It is known [12] that the following series holds for this Jacobi function 
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(16) 



being K(i) = 



1.3111028777 a constant. Then the mass spectrum of the theory in the limit of a very 



\) 4 Ho that we can recognize as the one of a harmonic oscillator. 



/0 \/l+sin 2 e 

large g is given by E n = (2n + 1) v 2) 

So, the main physical motivation to study our approach through Green functions for nonlinear equations is to have 
a deeper understanding in quantum field theories but the method is rather general and could find applications in a 
lot of other fields. 



IV. GREEN FUNCTION METHOD FOR NONLINEAR DIFFERENTIAL EQUATIONS 



In order to make our approach as clearer as possible, we consider the trivial problem of a Riccati equation 

y + y 2 = l (17) 
0. The solution is given by y(t) 



with the initial condition y(0) 
this equation being given by 



G(t) = 9(t) 



1 



1 + t 



tanh(i). A Green function is easy to compute for 

(18) 



So, let us consider the following small time expansion as a solution of the above Riccati equation 
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(19) 
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being a, b and c constants to be computed. In order to compute these constants we consider the equation we started 
with and compute all the derivatives till the order we are interested in. Then, we compare these derivatives with 
the one obtained through equation (|19|) fixing in this way the values of the constants to make them equal. So, from 
eq. lfTT))) one gets 



f 5 



y(t) w ln(< + 1) + a[- ln(t + l)+t]+ b[ln(t + 1) + — _*]+ c [- \ n (t + 1) + - + t - —] + 



(20) 



2 v ' 3 2 

and from this we can compute y(0), 2/(0), y(0) and so on. From the Riccati equation we have j/(0) = 0, y(0) = 1, 
?y(0) = and so on giving finally a = 1, b = — 1 and c = — 1 for our case yielding y(t) — t — y that are the first 
two terms of the Taylor series of the tanh(t), the exact solution of the equation. From this exercise we learn that the 
series p^|) is a small time series solution of the original equation and that the convergence may be really slow. It is a 
rather interesting aspect of this approach that the Green function method has such a way to be applied to nonlinear 
equations. The case we considered here is a rather trivial one but things are made more interesting for the case of a 
A</> 4 when we go to a numerical comparison. 

We want to apply the above approach to the case of equation (fT3"|) . So, let us seek a solution in the form (properly 
normalized by /Jo) 



sn 



- sn 
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(21) 
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FIG. 1: Comparison for a driving source j(t) = sm(27ri). 
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FIG. 2: Comparison for a driving source j(t) — exp(— t). 



and after computing derivatives of this equation and eq. pUl) we get easily a — ^ and b = — ^[j(0)] 2 . This gives 
the result we aimed for. We have got the proper expansion by Green function method of a solution to a nonlinear 
differential equation. What we want to see is how good is this approximation when compared to numerical results. 
This will be shown in the following section. 



V. NUMERICAL RESULTS 



In order to verify the quality of our approximation we solve the equation Q13[) for two different driving sources and 
take the coupling constant g = 1. Firstly we considered j(t) = sin(27rt) and the results are given in fig[T] In this case 
we can only have a first order correction as j(0) = 0. The agreement is very satisfactory till the end of the integration 
interval. 

The second case we considered j(t) = e _t permits to introduce another correction term but we notice no significant 
improvement due to the slow convergence of the approximation as can be seen from figdl 

The quality of the approximation depends on j(t) that can make very demanding the need for higher order correc- 
tions. 
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FIG. 3: Comparison between our approach and functional iteration method. 



VI. COMPARISON WITH OTHER METHODS 

There exist different techniques to manage non- linear differential equations. In order to have a proper comparison 
we have to limit our analysis to small time methods. In this sense, the most similar approach to ours is the functional 
iteration method (l3l | . This method proves to have a rapid convergence to a good approximant of the true solution 
when the equation is not stiff. This is exactly our case. So, let us consider the equation 



4>(t) + <f) 3 {t) = sin(27rf) 



(22) 



where dot means a time derivative. Functional iteration method implies that we solve the above equation iteratively. 
We assume 0(0) = and 0(0) = 0. We take as zero order iterate 4>o(t) = 0(0) = and then for the successive iterates 
we have 



$u+i(t) = -03(i)+sin(27rf) 



(23) 



starting with v — 0. Already at the second iterate we get a very good approximation to the true solution in the range 
we are considering. Then, we can compare this approximation with our method considering two terms. The results 
are presented in figEJ 

This result shows that functional iteration method has a faster convergence and at least another term should be 
computed with our approach to reach an identical precision in the required range. This result should also be expected 
on the ground of efficiency of iterative methods with respect to series solutions. So, in order to decide the proper 
method to use one should properly analyze the problem at hand. 



VII. CONCLUSIONS 



We have shown an approximation method to solve nonlinear differential equations using Green function methods. 
This method proves to be a small time expansion and the convergence in some cases may turn out really slow. The 
main point to be emphasized is the unexpected utility of this approach generally assumed to hold only for linear 
differential equations. This implies that a gradient expansion for nonlinear PDE can also be applied successfully 
and a quantum field theory obtained. In this latter case one should consider that a gradient expansion is a strong 
coupling expansion and then, information in this regime of the corresponding quantum field theory is given. This 
yields another method to approach these problems generally very difficulty to manage with analytical methods. 
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